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Abstract 

We consider a Lattice Gas model in which the sites interact via 
infinite-ranged random couplings independently distributed with a 
Gaussian probability density. This is the Lattice Gas analogue of 
the well known Sherrington-Kirkpatrick Ising Spin Glass. We present 
results of replica approach in the Replica Symmetric approximation. 
Even with zero-mean of the couplings a line of first order liquid-gas 
transitions occurs. Replica Symmetry Breaking should give up to a 
glassy transition inside the liquid phase. 
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It is well known [I]] that Ising Model is equivalent to Lattice Gas, a system 
denned in terms of occupation variables r taking values and 1. The Lattice 
Gas effective Hamiltonian is formally identical to Ising onef|TJ. A simple 
change of variables (a = 2r — 1) maps each of the two Hamiltonians into 
the other, provided that Ising external field is related to lattice gas chemical 
potential by h — J = ^/z, where J is Ising spin-coupling related to lattice 
gas site-coupling $ by J = This reflects into a simple relation between 
Ising free energy density and Lattice Gas pressure: p = h — | J — / . The 
two systems have therefore the same phase diagram and the same critical 
behaviour (real gases and Ising magnets are in the same universality class). 

For random systems 0] this whole argument breaks down because the 
relation between chemical potential and magnetic field involves the quenched 
couplings. As we shall see in the following this reflects in new and unexpected 
features for the phase diagram of the system. 

For Neural Networks the inequivalence between spin(±l) and occupation 
(0, 1) variables was already been pointed out and analyzed, see for example 
|H and references therein. 

Recently much effort has been devoted to develop a description of the 
structural glass transition ^, |J within the framework of disordered systems. 
All these models were however based on Ising Spin variables instead of Lattice 
Gas ones which would be more appropriate for a condensed matter system. 
For disordered systems the two kind of variables are not equivalent. To have 
a comparison term it would be useful to analyze the properties of a mean 
field disordered lattice gas model. 

We consider a system of N sites, an occupation variable is defined in 
each site k, can take the values or 1. The Hamiltonian of the system is 
taken to be formally identical to the SK's one||, ||. The interaction energy 
between two different (k and I) occupied sites is taken to be ipki and the 
system is coupled to some external source g. The total effective Hamiltonian 
is therefore: 

TV x 

h^t] = -9 22 T k- 9 famn . (l) 
k=i 1 i^k 

In magnetic language g would be the external field, while for a Lattice Gas g 
is the sum of the chemical potential, the kinetic contribution, and eventually 
an external force term. The infinite-ranged interaction energies {<pki} are 
taken to be quenched independent Gaussian variables with zero mean and 
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variance $/iV. In the following we shall take $ as our unit of energy and 
set $ = 1. We also set the Boltzmann constant equal to 1, as a consequence 
H, f, T, g and <fi are all dimensionless. Each <pki is taken to be equally 
distributed and therefore each site interacts with each other. 
For a given realization of the 0's the partition function is: 



to 



(2) 



We are interested, as usual when dealing with quenched disorder, to evaluate 
the averaged free energy density: 



mg) 



P[<j>\ InZ^dcf) 



(3) 



this will be done in the following using the replica approach^, 0, |J. 

We have to calculate the averaged n-th power of the partition function: 



z n = (z^y = J {z^) n p[(j>}d<i> . 

For integer n we get, after performing Gaussian integration: 



(4) 



Z n = E eX P 

to 



5>£tf + |^e(eW)' 

k a Ziv k<l \ a / 



we can reorder the exponent and obtain: 



Z n = E eX P 

to 



k \ a a 



a b 
k k 



Using Gaussian identities we rewrite (||) as: 



E ex P 

to 



Y[exp 

o,6 



N/3 2 V 4 V y 



(5) 



/?EUE^ a -7^E^ 



AN 



k 'k 



a.b 



X 



a.b 



'N/3'- 

47T 



2\ 2 



exp (- ^{3 2 Q 2 ab + \p 2 E Qab^rj/j dQ ab , (6) 
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reordering the exponentials and defining: 



T T 



In 



a.b 



we finally get: 



'N(3 2 

47T 



Je~ NA ^ d n2 Q 



The averaged free energy density is given by: 

f 



Af^oo n-*o pnJS 



In Z n (fi;g) 



(7) 
(8) 

(9) 
(10) 



In the thermodynamic limit the integral can be estimated by maximizing the 
integrand, this yields: 



fn(P;g) = Jim 



The extremum is determined from the saddle point equation: 

dA jP_ (3 2 J2t T a r b e~ HQ ^ 

Y Qab ~ T Ere-^M ~ ' 



dQ 



ab 



(11) 



(12) 



that may be rewritten as Q a b = {T a T b )Q. 

As a first stage we consider saddle points that are symmetric under the 
Replica Group 0. Setting Q ab = q + b5 a b we can write Y.abQab ra "T b = 
q (£„ r a f + b £ a r a and Eab Qlb = n(q + &) 2 + n(n - l)q\ substitution in (§) 
and extraction of the n — > limit then yields: 



/ = i /56(2g + 6)-(27r/3 2 )^ 



In 



X _|_ e 0(. a +*Vo) 



e a* 2 



(13) 



In equation ( |T3"D we set a = g + and the matrix elements satisfy the 
coupled equations: 



p = q + b = (2tt) 



1 + e 



e 2 cfe 
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q = (2tt)-3 / [l + e ^+ z Vz) 



-2 



dz . 



(14) 



As can be seen following the line of ||, || the physical significance of p and q 
is: 

p = —) q = W, (15) 



q = (r) 2 , 

a bar denotes the average over quenched 



where, following the notations of 
disorder. 

For 8 = we have p = \ and q — |, as we expect from their physical 
significance. In the high temperature regime we can solve ([14]) by expansion 
in powers of /3, this yields: 
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p = I + 1^ + I/? 2 

- + ~/3c H + g 2 

4 4 Py 16 U y 
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We have numerically solved equations (|14j) for several values of T and g. 
We find a line of first-order phase transitions. Such line ends in a second- 
order transition point, for T ~ 0.22 and g ~ —0.7, where the linear response 
function \ — dp/dg (the "susceptibility") is found to be divergent. 

Let us define 

7o = lim B(p — q) = lim 8b , 
for T — >• we have p = q + 7T and thus a = g + |7q. For /3 > 1 the 



solution of (|14 ) can behave in two different ways depending on which side of 
the first order transition we consider. We call go = —0.63633 the transition 
point at zero temperature. In the range g < go, qo (= q(T] <7)|t=o) an d 7o(fl0 
are identically zero (details will be presented elsewhere [|TT|| ), we find that all 
their temperature derivatives vanish for T — > 0, q vanishes as e^ 9 (g < 0) 
and 7 (= 8(p — q)) as /Je^ 9 . Elsewhere if g > go, they depend linearly on T, 
indeed we get: 



q = qo 



<?o7o(<?o + |7oao) 



(qo + ho®o) 2 + z7o(<?o - Oil] 



7 = 7o + 



ho(qo 



an 



(qo + ho^o) 2 + ho{qo 



org) 



T 



T 
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p = q + 7T . 

Next we look to thermodynamical functions, the internal energy and en- 
tropy densities are given by: 

u = -{g + lo)q 

s = ■ (16) 

The entropy fllDD is negative in the range g > go where 70 is different from 
zero, so we should expect the Replica Symmetry to be broken in this region. 
We stress that differently from SK case we have a region in which the Replica 
Symmetric Solution remains physical down to zero temperature. This can 
be related to some exact results on the ground state energy of the model [O]. 
The maximum absolute value of the zero temperature entropy is at g = g$ 
where it takes the value 0.101 and it strongly decrease for higher values of g 
(e.g. at g — 1 s = —0.011). This solution is clearly unphysical (for g > g ) 
at low temperature. This should be considered a clear signal @ of Replica 
Symmetry Breaking (RSB). 

Now we briefly compare our results with the ones of Sherrington and 
Kirkpatrick. In SK case the zero-field Hamiltonian has a global Z 2 symmetry 
and all thermodynamical functions are either even or odd in the external 
field. SK model can have a (field-driven) first-order phase transition only if 
a strong enough ferromagnetic part (i.e. a non-zero mean) is added to the 
random coupling. Such a transition is indeed related to the breaking of the 
global Z2 symmetry induced (as in the homogeneous case) by a ferromagnetic 
coupling. Moreover, at each value of external field, SK's Replica Symmetric 
Solution always becomes unphysical for low enough temperature. 

It is apparent how our picture is different from the usual one. We have 
considered the case of purely (zero-mean) random interactions and our Hamil- 
tonian has no Z 2 symmetry. The Replica Symmetric Solution of our system 
exhibits a line of first-order phase transition points ending with a second- 
order transition. This feature is expected to be robust to Replica Symmetry 
Breaking, as should be confirmed by the stability analysis [|l(], |ll|]. Two phases 
co-exist along the transition line and in one of them the Replica Symmetry 
should remain exact down to zero temperature. Following the conventional 
Lattice Gas interpretation this two phases should be called 'gas' and 'liquid'. 
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The gas-like phase is that which is replica-stable at all temperature. When 
Replica Symmetry Breaking sets up the liquid-like phase can give up to a 
glassy state. This idea will be put forward elsewhere j|TT| . 

I am very grateful to Prof. Giorgio Parisi for the help he gave in this 
work with his suggestions and experience, and also for his careful and critical 
reading of the original manuscript. I acknowledge Prof. Enzo Marinari and 
Dr. Felix Ritort for useful and stimulating discussion. I wish acknowledge 
INFN, "Sezione di Roma Tor Vergata" , for support received during this work. 
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